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ABSTRACT 

In this paper we present ground-based high-resolution spectroscopy of 51 Pegasi using 
CRIRES at the Very Large Telescope. The system was observed for 3x5 hours at 2.3 fim at 
a spectral resolution of R = 100,000, targeting potential signatures from carbon monoxide, water 
vapour and methane in the planet's dayside spectrum. In the first 2x5 hours of data, we find a 
combined signal from carbon monoxide and water in absorption at a formal 5.9a confidence level, 
indicating a non-inverted atmosphere. We derive a planet mass of Mp — (0.46 ± 0.02)Mjup and 
an orbital inclination i between 79.6° and 82.2°, with the upper limit set by the non-detection of 
the planet transit in previous photometric monitoring. However, there is no trace of the signal 
in the final 5 hours of data. A statistical analysis indicates that the signal from the first two 
nights is robust, but we find no compelling explanation for its absence in the final night. The 
latter suffers from stronger noise residuals and greater instrumental instability than the first two 
nights, but these cannot fully account for the missing signal. It is possible that the integrated 
dayside emission from 51 Peg b is instead strongly affected by weather. However, more data are 
required before we can claim any time variability in the planet's atmosphere. 

Subject headings: Planets and satellites: fundamental parameters, Planets and satellites: atmospheres, 
Techniques: spectroscopic 



1. Introduction 

The first discovery of an exoplanet around a 
main-sequence star, t he G-dwarf 51 Peg asi, was 
announced in 1995 ( Mavor and Queloa ). The 
presence of the planet was inferred from the mea- 
surement of a periodic Doppler-shift in the stellar 
lines, which was interpreted as being due to the 
motion of the star around the center of mass of the 
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star-planet system. Since then, the radial veloc- 
ity technique has been revolutionary in exoplanet 
research, resulting in hundreds of discoveries. 

51 Pegasi b is unlike any planet in our own So- 
lar System: with a minimum mass of ^ 0.5 Afjup, 
it orbits at only ^ 0.05 AU from the star, well 
within the orbit of Mercury (0.47 AU). At the 
time of discovery, the generally accepted theories 
of planet formation and evolution could not ex- 
plain the presence and survival of a Jupiter-size 
planet at such short orbital distance. This resulted 
in an intense debate on whether the radial veloc- 
ity signal of 51 Peg really was caused by an unseen 
planet. 

The first concern was that, since the orbital in- 
clination i was unknown and therefore the mea- 
sured mass was only a lower limit, 51 Pegasi b 
could actually be a low-mass star or brown dwarf 
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seen face-on. This hypothesis was rejected because 
of the low probabihty of this geometrical configu- 
ration (2.5xl0~^ for the hydrogen-burning limit 
of 0.08 Mp)), the stellar p rojected rotational ve- 
locity ()Francois et al. I ll996l) . the absence of an X- 
ray signal from the system ( Pravdo et al.l Il996) . 
and eventua lly the discovery of other 51 Pegasi- 
like o bjects ( Marcv and Butler 1996 : Butler et al.l 
19971 ). At the same time. iLin et al.l (|l996l ) showed 
that 51 Pegasi b could have formed much fur- 
ther away from the parent star, and subsequently 
migrated inward through tidal inter action with 
the p rotoplanetary disk. Moreover, iRasio etliD 
(|1996[ ) showed that the timescale for orbital decay 
via tidal dissipation was possibly longer than the 
main-sequence lifetime of the star. 

Another concern was that stellar radial pul- 
sations or hot spots in the surface could mimic 
the Doppler signal induced by a short-period gi- 
ant planet. This hypothesis was also found to 
be unlikely, because the slow stellar rotation and 
low chromospheric activity are inconsistent with a 
high-amplitu de sign al with a period of a few days. 
Nevertheless, iGravl ^997) performed a bisector 
analysis on the spectral lines of 51 Peg, concluding 
that the RV signal was possibly caused by a non- 
radial, unknown mode of stellar oscillation, rather 
than by a planet. However, subsequent bisector 
analyses at higher spectral resolution showed no 
variations (Haitzes_et_ al. 1997, 1998). In the mean- 
time, adding confidence to the planet hypothe- 
sis, the first exoplancts in eccentric orbits were 
found, making it difficult to e xplain the RV sig- 
nals with stellar mechanisms (jMarcv and Butler 
19961 ). Eventually, transits of the hot- Jupiter 



HP 209458 b were obse rved (jHenrv et al.l l2000t 
Charbonneau et al]|2000f) . finally settling the de- 
bate on the planetary nature of radial velocity 
measurements and opening the era of exoplanet 
characterization. 

Although RVs are very successful at detecting 
planets, they offer little information about the 
planet itself, whereas transits reveal the planet 
radius and the system inclination, which, when 
combined with the radial velocity curve, deter- 
mine the planet mass, its mean density, and a 
constraint on the planet's internal structure. Fur- 
thermore, starlight can filter through the planet 
atmosphere during transit, showing an imprint 
of its atomic and/or molecular gases. Atmo- 



spheric absorption increases the planet effective 
radius and the transit depth at a given wave- 
length. Therefore, by measuring the planet ra- 
dius as a function of wavelength a transmission 
spectrum can be constructed. This technique 
has led to the iden t ificati o n of atomic sodium 
(ICharbonneau et all l2002t ISnellen et all l200i 
iRedfield et al ."2008V potassium fei ng et al.ll201l[ ). 
hvdroee n (fvid al-Madiar ct al. 200^, carbon and 



water (iTinetti et al. 20071 


Swain et al.l 20081 


Desert et al.l 2009: Sine et al. 


20091 Gibson et al. 


2011). In addition, the planet is occulted bv the 



star around superior conjunction. Its light is tem- 
porarily blocked, allowing its flux to be measured 
directly by comparison with the out-of-eclipse to- 
tal system flux. Depending on wavelength, this 
reveals the planet's thermal em ission, possibly 



modulated by molecula r features (jCrillmair et al 



20081: ISwain et al.ll2009D, and/or reflected star light 



(jMazeh et al.ll2012l : lMislis and Hodgkinll2012l) . Fi- 
nally, continuous flux monitoring of a transit- 
ing sys tem can reveal the phase function of the 
planet (iKnutson et al.l[2009l : lcrossfield et al.ll20ld 



Knutson et al.ll2012[ ). 

Until recently, exoplanet atmospheric charac- 
terization was limited to transiting planets, leav- 
ing out most of the planets d iscovered via the RV 
technique. For the first time, Snellen et al. ( 2010h 
demonstrated that ground-based high-resolution 
spectroscopy and RVs can also be used for exo- 
planet characterization. They detected the ab- 
sorption signature of carbon monoxide in the 
transmission spectrum of HD 209458 b, using the 
CRIRES spectrograph on the ESO Very Large 
Telescope. Its spectral resolution of R~100,000 
allows the planet signature, which changes in 
Doppler shift during the observations due to the 
orbital motion of the planet, to be separated from 
the stationary telluric and stellar absorption lines. 
In addition to transmission spectroscopy, this 
method can be applied to dayside spectroscopy, 
probing the planet when it is almost fully illumi- 
nated, just before and after secondary eclipse. In 
this case the planet thermal emission, modulated 
by the atmospheric molecular signature, is directly 
detected and therefore does not require the planet 
to transit. This has led to the detection of car- 
bon monoxide absorption in the dayside spectrum 
of the non-transiting hot Jupiter r Bootis b, also 
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revealing its orbital inclination and planet mass 
(jBrogi et al.ll2012l : iRodler et al.ll2012l ). 



In this paper we apply the same method to the 
non-transiting planet 51 Pegasi b. The reminder 
of the text is organized as follows: in Section [2] 
we describe our observations and data reduction 
techniques. The resulting signal is described in 
Section|3l while Section |4] extensively discusses the 
non-detection in one of the three nights of obser- 
vations. Assuming that the signal in the first two 
nights is genuine, Section [5] presents the derived 
parameters of 51 Pegasi b and of its atmosphere. 
Finally, a recap of our main findings and the fu- 
ture perspectives of high-resolution ground-based 
spectroscopy are presented in Section IB] 

2. Observations and Data Reduction 
2.1. Telescope and instruments 

High-dispersion infrared spectra of 51 Pegasi 
(K = 3.91 mag) were taken with the ESO Very 
Large Telescope (VLT) as part of the Large Pro- 
gram 186.C-0289, aimed at the atmospheric study 
of the brightest transiting and non-transiting hot 
Jupiters visible from Cerro Paranal. The system 
was observed in 2010 on the nights of October 
16, 17 and 25 at low airmass, using the CRyo- 
genic Infra-Red Echelle Spectrograph (CRIRES, 
Kaufi et alll2004[ ). mounted at the Nysmith-A fo- 
cus of the VLT Antu. The spectra are imaged on 
four 1024x512 pixels Aladdin II detectors, sepa- 
rated by small gaps of about 100 pixels. At a 
resolution of R ^ 100,000, 904 spectra were col- 
lected during a total of 16 hours of observations 
in the wavelength range 2.287-2.345 /im, which 
covers the 2-0 R-branch of carbon monoxide. The 
Multi- Applicatio n Curvature Adaptive Optic sys- 
tem (MACAO, lArsenault erall l2003l ) was em- 
ployed to maximize the throughput of the 0.2" slit. 
During each night, the target was observed with- 
out interruption while nodding along the slit by 
10" between consecutive spectra, according to an 
ABBA pattern, to allow a proper background sub- 
traction. A standard set of calibration frames was 
taken in the morning after each night of observa- 
tion. 



2.2. Basic data analysis 

The initial data reduction was performed us- 
ing the CRIRES pipehne LlLcQ. Each set of 
AB or BA spectra was flat-fielded, corrected for 
bad-pixels and non-linearity and combined in or- 
der to subtract the background. One-dimensional 
spectra were subsequently e xtracted with the op- 
timal extraction algorithm ( Hornd 1986 ). result- 
ing in 452 spectra. The subsequent data analy- 
sis was performed with purpose-built IDL scripts. 
For each night and each detector, the spectra were 
handled as a matrix with wavelength (or pixel 
number) on the a;-axis and time (or phase) on 
the y-axis, and treated separate ly, f ollowing the 
proced ure in [Snellen et al. I (|2010l) andlBrogL et al.l 
(|2012[ ). This is necessary because each night of 



observations has different atmospheric and instru- 
mental conditions, and because each detector is 
read out using a different amplifier, with particu- 
lar characteristics that need to be modeled inde- 
pendently. 

In each spectrum, bad pixels and bad regions 
were corrected through spline and linear inter- 
polation respectively. The identification of bad 
pixel/regions was first done by eye on each ma- 
trix of data. Residual outliers were removed itera- 
tively during the whole reduction sequence. Sub- 
sequently, each spectrum was aligned to a com- 
mon wavelength scale. The difference between 
the centroids of the telluric lines and their aver- 
age value across the entire series was fitted with 
a linear function in pixel position, and the result- 
ing shift was applied via spline-interpolation. Af- 
ter the alignment the typical scatter in the resid- 
ual position of the line centroids is on average less 
than 0.1 pixels. The common wavelength solution 
was determined by fitting a second-order polyno- 
mial to the pixel positions of the lines in the av- 
erage spectrum, as function of their correspond- 
i ng wavelengths listed in the HITRAN database 
( Rothman et al. I I2OO9I) . Finally, each spectrum of 
the series was normalized by its median value. 

2.3. Stellar subtraction 

Compared to our previous analysis of HD 
209458 and r Bootis, an additional step - the 
removal of the lines in the stellar spectrum - was 
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required. This is because 51 Pegasi is a compara- 
tively cooler star with an effective t emperature of 
Tes = (5790 ± 44) K (|Fischer and V alcnti 2005), 
and thus several spectral lines, in particular those 
of CO, are also present in the stellar spectra. Due 
to the changing barycentric velocity of the obser- 
vatory during the night (Awobs ^0.5-1 km s~^), 
mainly caused by the rotation of the Earth, the 
stellar lines shift by 0.3 - 0.6 pixels in wavelength 
during our observations. If the stellar spectrum 
is not subtracted from the data, the procedure 
used to remove the telluric contamination (see 
Section [2^ would produce strong residuals at the 
position of the stellar lines. Since the effective 
temperature of 51 Pegasi is very similar to that 
of the Sun, a stellar template was created from 



a hig h-resolution solar spectrum (jAbrams et al 
19961 ). The position, amplitude and width of the 



solar lines were fitted with Gaussian or Lorentzian 
profiles (where appropriate) and a model spec- 
trum was built from the list of the fitted lines and 
convolved to the CRIRES resolution. Each frame 
was then divided through the normalized template 
spectrum, shifted in wavelength according to the 
radial velocity of the star 51 Pegasi, and multi- 
plied by a factor of 0.95 in order to match the line 
depths relative to those of the Sun. The need for 
a scaling factor is possibly due to differences in 
spectral type, metallicity and/or rotational veloc- 
ity between the two stars. 

2.4. Removal of the telluric contamination 

The most dominant feature in the spectra is 
telluric absorption, which has a fixed pattern in 
wavelength but is variable in depth due to changes 
in airmass, atmospheric conditions and instrumen- 
tal response. In order to minimize the telluric sig- 
nature, the logarithm of the fiux in each column 
of the spectral series was first linearly fitted with 
geometric airmass and divided out. This greatly 
reduced the amplitude of the telluric features, 
although residuals at the 1-5% level were still 
present at the positions of some of the strongest 
absorption lines, due to second-order effects. For 
example, possible changes in the water vapour 
content during the night will generally not follow 
airmass. The behaviour of the residuals as func- 
tion of time was therefore measured in a few strong 
lines and the flux in each column of the matrix 
was fitted by a linear combination of those mea- 



surements, computed via linear regression. In ad- 
dition, residual large-scale gradients in the matrix 
were removed by applying a high-pass filter to each 
spectrum of the series. Finally, each column was 
divided by the square of its standard deviation, 
in order to scale down the noisy parts of the spec- 
tra according to their signal-to-noise. These would 
otherwise dominate the cross-correlation analysis 
performed later. 

2.5. Atmospheric models 

In order to be able to combine the spectral lines 
of 51 Pegasi b via a cross-correlation technique, 
model spectra of the planet atmosphere were com- 
puted assuming an average dayside vertical pro- 
file in hydrostatic equilibrium. The atmosphere 
is assumed to be isothermal for pressures p > pi 
and p < P2, at temperatures Ti and T2 respec- 
tively. In between the two pressure levels, a con- 
stant lapse rate [d\og{p)]/dT is assumed. Inverted 
and non-inverted T /p profiles are modelled by as- 
suming T2 > Ti and T2 < Ti respectively. Such 
a vertical structure is the simplest possible for a 
planet atmosphere, and although it is probably 
not representative for the entire vertical extent, it 
is a reasonable assumption in the limited range 
of temperature and pressure probed by our obser- 
vations. More importantly, given the degenera- 
cies between vertical profiles and molecular abun- 
dances (see Section [5^ . a more complex descrip- 
tion of the planet atmosphere would be beyond the 
purposes of this study. Molecular species were as- 
sumed to be vertically mixed, with the dominant 
sources of opacities at 2.3 micron being carbon 
monoxide, water and possibly methane. Opacities 
for CO a nd H2O were taken fr om the HITEMP 
database ( Rothman et al. 2010l ). while those for 



CH4 were taken from HITRAN (jRothman et al 
2OO9I ). A grid of models was computed assum- 



ing volume mixing ratios (VMRs) between —4.5 < 
logio(VMR) < -3.0 in steps of Alog^o (VMR) = 
0.5 for the three molecular species, at pressures of 
pi = (1, 0.1, 0.01) bar and p2 = (10-^ lO"*, 10"^) 
bar, and at three base temperatures of Ti = 
(1000, 1250, 1500) K. These are representative of 
planet equilibrium temperatures in various cases 
of energy redistribution. T2 was fixed to 500 K 
for non-inverted profiles, and to 1500 K for in- 
verted profiles. All possible combinations of pa- 
rameters were employed. The radiative transfer 
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Fig. 1. — Example model spectra for the atmosphere of the planet 51 Pegasi b, showing the effects of different 
base temperatures Ti (left panel) and different Volume Mixing Ratios of water and carbon monoxide (right 
panel). Details about these models are given in Section [23] 



in the planet atmosphere was computed in the 
plane-parallel approximation and typical resulting 
models are shown in Figure [TJ As explained in 
Section 15. 2[ in order to set lower limits on the 
VMRs of molecules, we also assumed adiabatic 
lapse rates. In this case, p2 is set by the altitude 
at which the atmosphere reaches T2 = 500 K. 

2.6. Cross correlation 

In our observed wavelength range, the expected 
spectral signature of 51 Pegasi b consists of tens of 
carbon monoxide lines, plus a possible additional 
contribution from water vapor and/or methane, 
with a typical depth per-line of few parts in 10^ 
with respect to the stellar continuum. Conse- 
quently, the individual lines will be in any case 
deeply buried in the noise, as each spectrum has a 
signal-to-noise ratio of typically 500. In order to 
detect the planet signature, the signals from all the 
lines a re combined with a cross-correlation tech- 
nique (jSnellen et al.ll20ld iBrogi et al.ll2012[ ). In- 
stead of directly cross-correlating the atmospheric 
models with the data, a template spectrum is con- 
structed from the model by subtracting the base- 
line, clipping out the weakest lines and convolv- 
ing the result to the CRIRES resolution. This 
is done because any broad-band signal has been 



removed by our data reduction technique. In ad- 
dition, the line profiles in the model should match 
the instrumental profile in order to maximize the 
cross-correlation signal. Each spectrum is cross- 
correlated with the template using a lag-vector 
corresponding to radial velocities between -250 
and -1-250 km s~^, in steps of 1.5 km s~^, and the 
result for each detector and each night is stored 
in a matrix with x- and y-dimension representing 
radial velocity and frame number (phase) respec- 
tively. Subsequently, the matrices of the four de- 
tectors are summed with equal weighting, result- 
ing in a single cross-correlation matrix per night 
of observation. As an exception, detector 4 was 
discarded from the analysis during the third night 
(October 25, 2010), as explained in Section [4.11 
For models with non-inverted (inverted) T /p pro- 
files, a positive cross-correlation signal indicates 
absorption (emission) in the planet atmosphere. 
At this stage the signal-to-noise is insufficient to 
detect the peak of the cross-correlation function 
per spectrum so we cannot derive the planet ra- 
dial velocity curve on a frame-by-frame base. In- 
stead, the cross-correlation signal needs to be com- 
bined over time by considering a range of maxi- 
mum planet radial velocities {K^) corresponding 
to the possible system inclinations. In all cases 
we assume a circular orbit. For each value of -ftTp, 
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Fig. 2. — Left column: Total cross-cross correlation signal for each night of observations, without removing 
the stellar absorption lines. Right column: Same as the left column, but after removing the stellar spectral 
features. In all panels a linear color scheme running from -7(t (black) to +7a (white) is adopted. Crosses 
indicate the planet position measured by combining the nights of October 16 and 17. Black-dashed lines 
encompass regions of the parameter space significantly affected by stellar residuals. These result in extra 
cross-correlation noise potentially masking the planet signal in the third night, as discussed in Section [4. II 
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the planet radial velocity curve was computed and 
utilized to shift each cross-correlation function to 
the rest frame of the planet, via linear interpo- 
lation. Finally, the shifted cross-correlation series 
was summed over time, resulting in the total cross- 
correlation signal as function of systemic radial ve- 
locity (Ksys) and maximum planet radial velocity. 
If 51 Pegasi b is detected, we should find a peak 
in the cross-correlation at the known systemic ve- 
locity of the system (14ys = —33.25 km s~^), for 
the best matching inclination. Any spectral signa- 
ture which is not accelerating with respect to the 
observer will appear at Kp — 0, including resid- 
ual telluric and stellar absorption features. More 
importantly, since the cross-correlation signal is 
summed in time by assuming a range of inclina- 
tions, strong telluric and stellar signals will also 
show residuals at Kp ^ 0. These will be shifted in 
systemic velocity in a non-trivial way, depending 
on the phase range of the observations (see e.g. 
the left column of Figure ^ . This is why telluric 
removal and stellar subtraction are crucial steps of 
our data analysis. 

3. The cross-correlation signal 

The total cross-correlation signal, computed as 
explained in Section 12.61 is shown in Figure [2] 
for each of the three nights of observations, us- 
ing the planet model spectrum that best matches 
the data. This is selected from our large grid of 
models as explained in Section r5.2l and corresponds 
to atmospheric absorption (non-inverted T/p pro- 
file) of carbon monoxide and water, with VMRs of 
10""* and 3 x 10^'' respectively. On the nights of 
October 16 and 17, 2010, candidate signals from 
51 Pegasi b are detected at the orbital velocities 
and systemic velocities listed in Tabled] Residuals 
from stellar lines appear as extra cross-correlation 
noise in a region of the parameter space approx- 
imately enclosed within the black dashed lines in 
Figure [2] The night of October 25 does not show 
any significant planet signal, and we discuss pos- 
sible reasons for this in Section |4l 

Figure [3] shows the cross-correlation signal of 
the nights of October 16 and 17 combined. A 
peak in the total cross-correlation signal is de- 
tected at the orbital radial velocity of Kp — (134.1 
± 1.8) km s"^. Its significance of 5.9a is com- 
puted by dividing the peak cross-correlation value 



by the standard deviation of the noise. This 
points to a high orbital inclination, as discussed 
in Section 15.11 We discuss how this signal com- 
pares to that obtained using models with pure 
water and pure carbon monoxide in Section 15.21 
A phase shift of Aip — 0.0095 needs to be ap- 
pli ed to the orbi t al par ameters of the system given 
bv lButler et"aI1 (|2006l ) in order to match the sig- 
nal to the known systemic velocity of 51 Pegasi 
(i^sys = —33.25 km s~^). This shift is well within 
the la uncertainty range of A(p = ±0.012 derived 
from the same orbital solution. 
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Fig. 3. — Total cross-correlation signal obtained 
by combining the nights of October 16 and 17 and 
by applying a global phase shift of Aip = 0.0095. 
The planet signal is recovered at the systemic ve- 
locity of 51 Pegasi {Vsys = —33.25 km s~^, white 
dashed line) and at a maximum planet radial ve- 
locity of Kp = (134.1 ± 1.8) km s~^ Note that 
there is an uncertainty in linking the measured Kp 
(left vertical axis) to the orbital inclination i (right 
vertical axis). This is due to the uncertainties in 
the stellar mass (see Section [5TT|) . 



Since one third of the data do not show the 
planet signal, it is particularly important to ro- 
bustly assess the statistical significance of the sig- 
nal in the first two thirds of the data. For this pur- 
pose, the properties of the cross-correlation noise 
were studied by comparing the distribution of the 
cross-correlation values with that of a Gaussian 
distribution. This is shown by the left panel of 
Fig. 131 In addition, in the right panel the quantiles 
of the observed distribution are plotted against 
the quantiles of a Normal distribution (Q-Q plot, 
Wilk and Gnanad csikan 1968). In both cases, the 
two distributions agree within ^4(T, a limit set by 
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Fig. 4. — Left panel: the distribution of the cross-correlation values for the nights of October 16 and 17, 2010 
is shown in logarithmic scale, with la error bars corresponding to the square root of the bin occurrences. 
Overplotted in grey is the expected distribution for Gaussian noise. Right panel: the quantiles of a Normal 
distribution are plotted against the quantiles of our measured cross-correlation values (Q-Q plot). In both 
panels, the observed and theoretical distributions agree across the entire range which can be probed (±4cr), 
which is set by the number of cross-correlation points available. 



the finite number of points available. It is not 
surprising that the cross correlation has Gaussian 
noise properties, as non-Gaussian noise originating 
from a variable signal-to-noise is avoided by nor- 
malizing the data by the square of their standard 
deviation (see Section 12. 4|) . Moreover, system- 
atic noise is also strongly scaled down by the fact 
that the template spectrum contains tens of well- 
spaced molecular lines, and therefore the cross- 
correlation spans large portions of the data. 

In Figure [5] we compare the distribution of the 
cross-correlation values inside the planet trail (i.e. 
those around the planet radial velocity for the 
derived Kp) with that of the points outside the 
planet trail. The in-trail distribution is shifted 
towards higher values than the out-of-trail distri- 
bution. In order to quantify thi s, a Welch t-t est on 
the two samples was performed ( Welchl[l947 ) . with 
the null hypothesis that they are drawn from the 
same parent distribution. Based on the results of 
the test, the null hypothesis is rejected at the 5Aa 
confidence level, in line with our initial estimate 
based on pure signal-to-noise calculation. More- 
over, Fig. ini shows the Welch t-test performed for 
the full range of T^ys and Kp of Fig. [31 Remark- 
ably, no other region of the parameter space shows 



significant cross-correlation signal, except for two 
slanted lines converging towards the planet signal. 
This is expected from a real planet signal, because 
the method used for summing the cross-correlation 
functions in time naturally produces signals fad- 
ing away from the planet position. For this reason 
two or more observations taken at different phases 
are much more effective in constraining the planet 
parameters than a single observation. 

4. The non-detection on the third night 

As explained in Section [3l the third night of 
observations does not show any sign of an ab- 
sorption signal from the dayside of 51 Pegasi b. 
This is a concern because of the statistical robust- 
ness of the detection during the first two nights 
(see Section [3|), which makes it unlikely that time- 
correlated noise could mimic the planet signal. 
Therefore, it is important to investigate the possi- 
ble causes of this discrepancy. We discuss below i) 
possible issues with the instrument, data and/or 
data analysis for the third night, and ii) possible 
astrophysical causes. 
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Fig. 5. — Normalized distribution of the cross- 
correlation values out of the planet trail (black 
line) compared to the same distribution for the 
in-trail cross-correlation values. The in-trail dis- 
tribution is centered to higher values than the 
out-of-trail distribution at a 5Aa confidence level 
(see Figure [5]). The error bars corresponds to the 
square-root of the bin occurrences. 



4.1. Issues with instrument, data and data 
analysis 

The night of October 25 experienced a rela- 
tively good seeing and a longer coherence time 
than the other two nights of observations, result- 
ing in better performances of the MACAO system 
and therefore in a higher signal to noise compared 
to the first two nights (by 5-15% per spectrum). 
However, the stability of the CRIRES spectro- 
graph was significantly worse, causing the spectra 
to drift by ~1.5 pixels in the dispersion direction, 
whereas during the other two nights these shifts 
were only 0.2-0.5 pixels. The cause of this drift is 
likely associated with thermal instabilities of the 
cryogenic system. The temperature of the pre- 
optics system, grating stabilizer and grating itself 
changed by about 1.5 K during the 5 hours of ob- 
servation on October 25, whereas during the other 
two nights it was stable at the < 0.05 K level. This 
drift in temperature is indeed correlated with the 
drift in spectral position on the detector. Since the 
spectra need to be aligned to a common reference 
frame (see Section 12. 2|) , the wavelength drift has 
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Fig. 6. — Results of a generalized t-test on the 
null hypothesis that the in- and out-of-trail cross- 
correlation values are drawn from the same parent 
distribution. This test is performed for the same 
combination of systemic velocities and planet or- 
bital velocities as in Figure[3] It indicates that the 
null hypothesis Hq can be rejected at the 5.4cr con- 
fidence level for the planet signal. Note that there 
is an uncertainty in linking the measured Kp (left 
vertical axis) to the orbital inclination i (right ver- 
tical axis). This is due to the uncertainties in the 
stellar mass (see Section [5?T|) . 



a negative impact on the quality of the data. We 
were forced to discard detector 4 from the analy- 
sis of the third night, because it suffers from gain 
variations between neighbouring columns (the so- 
called odd-even effect). While with sub-pixel shifts 
this is not a significant problem (spectral features 
roughly fall on the same pixel), for shifts of one 
pixel or more the odd-even effect introduces strong 
time-correlated noise in the data of up to several 
percent. In addition, general flat- field inaccuracies 
are affecting the data more significantly. 

Despite the instrumental issues during the third 
night, the expected significance of the detection in 
the third night is ~ 3.8ct, based on the measured 
signal-to-noise, the total number of spectra, and 
the fraction of data lost because of the exclusion 
of detector 4. This means that, taking into ac- 
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count statistical fluctuations, there is a probabil- 
ity of 20%, 10% and 5% to measure a signal less 
than 2.5a, 2.2a and 1.8ct respectively during the 
third night. 

Another potential issue concerns residual stel- 
lar spectral features. As explained in Section 12.61 
and shown in Figure [21 even after subtracting a 
model stellar spectrum from the data, residual 
cross-correlation noise is seen. This is likely due to 
changes in the spectral response profile during the 
night and to inaccuracies in the alignment. This 
is particularly an issue during the third night, be- 
cause the system was observed at phase (p ~ 0.5, 
when the planet spectral lines overlap with those 
of the star (i.e. planet and star have almost the 
same radial velocity). In contrast, during the 
nights of October 16 and 17 the planet and the 
star have a different radial velocity, causing this 
residual cross-correlation noise to be located in re- 
gions of the parameter space far away from the 
planet position, and therefore not significantly in- 
fluencing the measured signal. The effects of the 
stellar residuals are shown in Figure [21 before and 
after the subtraction of the stellar spectrum (left 
and right columns respectively). During the third 
night of observations the expected planet signal 
falls clearly on top of the residual stellar cross- 
correlation (bottom- left panel) . After the subtrac- 
tion of the stellar spectrum residual signal is still 
visible (bottom-right panel). It is possible to filter 
out these stellar residuals by masking 5-7 pixels 
in the cross-correlation matrix, around the stellar 
radial velocity. In this way another ^ 40% of the 
planet signal is removed, as estimated by injecting 
a model spectrum in the data at 3x the nominal 
level and at the previously determined Kp and re- 
running the data analysis. Even after masking, 
no significant planet signal is detected. Although 
stellar residuals cannot explain the complete ab- 
sence of a planet signal, they certainly damp part 
of it. This is another important factor that differ- 
entiates the third night from the other two. 

4.2. Astrophysical causes 

We also explored which possible astrophysical 
effects may cause the non-detection during the 
third night. First of all, the planet is observed 
around phase ip = 0.5, meaning that if it was 
a transiting system the planet would be eclipsed 
by the star, potentially explaining why no signal 



was detected. However, photometric monitoring 
of the 51 P eg system sh ows that the planet doe s 
not transit (|Henrv et al.l[l9 97: Wal ker et al.ll2006[ ). 
For a system to exhibit secondary eclipses without 
transits, we calculated that the orbital eccentricity 
needs to be at least e = 0.13 for a stellar radius of 
i?s = (1.237±0.047)i?Q (|van Belle and von Brau3 



2OO9I ) and a planet radius of 1 i?jup- This is clearly 



excluded by the ra dial velocity data of the star 



( Butler et al.ll2006[ ) 



Small deviations from the best-fit orbital solu- 
tion to these radial velocity data can infiuence the 
radial velocity and timing of the expected planet 
signal, which could cause the non-detection during 
the third night. All the orbital solutions available 
in the literature agree over the fact that the ra- 
dial velocity data of 51 Pegasi are equally well fit 
by a circular orbit as by a slightly eccentric orbi t 



by a circular orbit as by a sligntly eccentric orbi t 
(e.g.lMavor and Quelo3ll995t butler et al.ll2006[ ). 
iButler et all ( 20061) obtained a best solution with 
e = (0.013 ± 0.012) and w « 58°. We tested vari- 
ous eccentric solutions in order to retrieve the to- 
tal cross-correlation signal and at the same time 
the of the orbital solution to t he stellar radial 
velocity data (jButler et al.ll2006l ). Assuming an 



eccentricity of e = 0.025 and w ~ 180°, which is 
more than 2 a _away from the best-fit solution of 
Butler et al. ( 2006[ ). it was only possible to shift 
the signals of the third night by ^ 4 km s~^, 
whereas the closest cross-correlation peak in the 
third night can be aligned to the other two only if 
shifted by about 10 km s~^. This means it is not 
possible to explain the non-detection by assuming 
an eccentric solution. 

When discussing astrophysical causes, there is 
also the possibility that the signal from the at- 
mosphere itself is changing, either due to an ori- 
entation effect or due to to temporal changes. 
At orbital phase ip — 0.5, the orientation of the 
planet is almost precisely on the sub-stellar point, 
meaning that a relatively hotter atmosphere is ob- 
served. Although the presence of an inversion 
layer in hot- Jupiter atmospheres is still not well 
unde rstood, stellar irradiation is thought to play a 
role ( Fortnev et al. 12008,') . However stellar activ- 
ity ( Knutson et al. 120101) and diff erent elemental 
abundances ('M adhusudhanI 2012 ) have also been 



suggested as possible causes. Higher temperatures 
mean that some molecular species remain in the 
gas phase, and are capable of absorbing radia- 
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tion high up in the atmosphere, causing a ther- 
mal inversion and spectral features in emission. 
If the atmosphere of 51 Pegasi b would be par- 
tially inverted around the sub-stellar point, this 
may cancel out part of the signal. In the limit 
of a nearly isothermal dayside atmosphere, as it 
seei ns to be the case for some hot Jupiters (see 
e.g. lCrossfield et al.ll2ni2l) . no signal would be de- 
tected. However, the illumination of the visible 
planet disk changes only by ~20% during these 
observations, requiring a large temperature gra- 
dient as function of longitude and a significant 
contribution in absorption from the night-side, 
which is in contrast with an inefficient heat trans- 
port required by the presence of a much hotter 
sub-stellar spot. Moreover, such a large gradi- 
ent would also show up in our data as a phase- 
dependent cross-correlation signal, something that 
we do not detect. This leaves only the possibility 
of temporal variations, the details and timescale of 
which are largely unknown even for hot Jupiters. 
Although theoretical studies predict weather to 
be important in hot Jupiters and brown-dwarfs 
(jShowman et al.ll2Q12r ). with only three nights of 
observations we are reluctant to attribute the ab- 
sence of the planet signal to exoplanetary weather. 

5. Characterization of exoplanet 51 Peg b 

In the previous sections we showed that our ob- 
servations resulted in an apparent planet signal in 
the first two nights, but no detection in the third 
night. We have discussed several potential causes 
of the non-detection, including genuine astrophys- 
ical causes, none of which was robustly proven to 
be associated to the lack of signal. Even though we 
cannot confidentially exclude the possibility that 
the signal observed in the first two nights is due 
to a yet unknown and unfortunate instrumental 
effect, in this Section we assume that the signal is 
genuine and proceed with the characterization of 
51 Pegasi b. 

5.1. Planet mass and system inclination 

The known radial velocity amplitu de of the host 



of Ms = (1. 5 lb .04) Mq, as determined by 
Takeda et all (|2007[ ). it translates into a planet 
mass of Mp = (0.46 ± 0.02) Mjup. The orbital 
velocity of the planet can be computed through 
Kepler's Third Law, and it is equal to Vp = 
(133.7 ± 1.8) km s"^ The sine of the inclina- 
tion, sini = Kp/Vp = (1.003 ± 0.019), implies la, 
2a and 3a lower limits on the orbital inclination of 
79.6°, 74.7° and 71.0° respectively. Observations 
conducted with the MOST satellite show that the 



star, Ks = (55.94 ± 0.69) m s'^ (|Butler et al 



20061 ) ■ combined with our measured planet ra- 
dial velocity amplitude of Kp = (134.1 ± 1.8) km 
s~^, gives a star/planet mass ratio of Mg/Mp = 
Kp/Ks = (2397 ± 44). With a stellar mass 



planet is not transiting (^Walker et al.l 120061 ). Us- 
ing the planet and stellar radii given in Sect. 14.21 
this sets an upper limit on the system inclination 
of z„,ax= (82.2 ±0.3)°. 

5.2. The atmosphere of 51 Pegasi b 

The first information which is derived from the 
cross-correlation signal is that the atmosphere of 
51 Pegasi b is not thermally inverted in the pres- 
sure range probed by these observations. This is 
because significant planet signal is only observed 
when using models with non-inverted T/p pro- 
files (atmospheric absorption). In addition, sig- 
nificant cross-correlation signals are obtained not 
only from models with pure carbon monoxide, but 
also from models with pure water, at a < 5.7ct and 
< 4.1tT level respectively. Details of these single- 
molecule detections are given in Table [2] No sig- 
nificant signal from a pure methane model is ob- 
served. This suggests that, while in our previous 
study of r Bootis b the cross-correlation signal was 
entirely due to CO, here there is a significant con- 
tribution from H2O lines. In order to take this 
into account, another set of atmospheric models 
was computed by including both molecules, for the 
same range of pressure and temperature as in Sec- 
tion [231 Further information about the structure 
of the planet atmosphere can be derived by com- 
paring the cross-correlation signal retrieved by us- 
ing this additional grid of atmospheric templates. 
Note that, since the data analysis removes broad- 
band spectral features such as the slope of the 
planet continuum and the wings of broad lines, the 
cross correlation is only sensitive to narrow spec- 
tral features, such as the flux ratio between the 
planet continuum and the core of the molecular 
lines. The comparison between models and data 
can be done in the following two ways. Firstly, 
different model spectra, exhibiting varying line ra- 
tios, can be cross-correlated with the data. The 



11 



model spectrum that matches the data the best 
win show the the strongest signal. Secondly, the 
amplitude of the model spectrum can be compared 
with the dayside spectrum of 51 Pegasi b by sub- 
tracting the model spectrum from the data early 
on in the analysis chain, before telluric removal. 
This comparison must account for the dilution by 
the starlight, which results in the following nor- 
malized planet spectrum: 



.(A) = 



model 



(A) 



Rp 
Rs 



(1) 



where -Fmodci is the model spectrum, Fgtar the 
stellar spectrum (for which we use the NextGen 
low-resolution model for a 51 Peg-like star with 
Tgff = 5800 K, logg = 4.5, and solar metallic- 
ity) and Rp and Rg the planet and stellar radii 
respectively. This scaled model is reversed and 
injected in the data at the measured planet po- 
sition, before telluric removal. As a result, the 
cross-correlation signal will be cancelled out when 
the model matches the data. However, the fact 
that the planet is not transiting the host star adds 
significant uncertainty to this second test, because 
Rp is unknown. Where not specified, 51 Pegasi b 
is assumed to have the size of Jupiter. 

It is important to realize a degeneracy exists 
between molecular abundances and atmospheric 
lapse rate. The higher the VMR, the higher the 
altitude at which a molecule absorbs. This gener- 
ally translates to a larger temperature difference 
between the layer in which the planet continuum 
is formed and the altitude at which a molecule ab- 
sorbs in the core of a line. However, this difference 
is also related to the atmospheric lapse rate, i.e. 
the rate at which the temperature changes with 
pressure. Since the temperature difference deter- 
mines the contrast between the continuum and the 
core of the lines, the observed strength of the high- 
resolution planet signal is affected by this degener- 
acy, which cannot be solved with a single observa- 
tion in a limited wavelength range. In this work, 
the lapse rate was fixed to the adiabatic value (i.e. 
the maximum allowed in a planet atmosphere, cor- 
responding to the onset of convection), and the 
base temperature Ti was chosen to be roughly the 
planet equilibrium temperature (see Sect. 12.51 for 
details). This results in a maximum planet signal 
for a given molecular abundance and therefore sets 
lower limits in the molecular VMRs. 



Among the large set of models generated as ex- 
plained in Sect. 12. 5| we selected the best- matching 
dayside model spectrum of 51 Pegasi b by re- 
quiring, at the same time, the highest signal- 
to-noise of the total cross-correlation signal and 
zero residual cross-correlation signal after sub- 
tracting the model spectrum from the data. A 
range of ±lcr around these two values is explored 
in order to allow noise fluctuations. This cri- 
terium resulted in a model with Ti = 1250 K, 
P2 = 10-^ VMR(CO) = 10-^ VMR(H20) = 
3xl0~*, meaning that roughly 40% of the sig- 
nal is due to H2O molecular lines. This is not 
surprising, given the relatively lower equilibrium 
temperature of 51 Pegasi b as compared to t 
Bootis. Below 1500 K carbon monoxide starts 
to be significantl y less abunda nt than water for 
solar C/0 ratios (jMadhusudha n 2012). Our de- 
tection of water together with carbon monoxide, 
assuming the signal is genuine, is therefore in line 
with a relatively cooler hot Jupiter. Note how- 
ever that the relative abundance of water with 
respect to carbon monoxide is only weakly con- 
strained by these observations, with the la range 
being 0.05 < VMR(H20)/VMR(C0) < 140. 

Note that, intuitively, atmospheric models con- 
taining both water and carbon monoxide would 
produce a cross-correlation signal which is roughly 
the quadrature sum of that obtained from mod- 
els with pure H2O or pure CO, combining in this 
case to ^ 7(7 of significance. However, in prac- 
tice, this does not happen, because the two cross- 
correlation signals cannot be considered to be in- 
dependent. Superposition of CO and water lines 
(see e.g. Figure [U right panel), and a mismatch 
between the real planet spectrum and the model 
templates used, can cause the total signal to de- 
viate from a simple quadrature sum of that of the 
individual gases. Also, the combined model redis- 
tributes the noise content of the data. This means 
that the combined signal may deviate from the 
value expected from the individual gases with a 
probability distribution having a standard devia- 
tion of 1(7. 

6. Conclusions and future prospects 

By observing the orbital motion of 51 Pegasi 
b directly using high-resolution spectra, we have 
likely been able to disentangle the spectral finger- 
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prints of the planet's atmosphere from teUuric and 
stehar contamination, and reveal significant ab- 
sorption from carbon monoxide and water vapour 
in its dayside hemisphere using the first two nights 
of observations. This was achieved by studying the 
cross-correlation of the data with a model planet 
atmosphere. In this way the radial velocity of the 
planet is measured directly, allowing us to deter- 
mine the planet / star mass ratio and the system in- 
clination. This analysis suggests that 51 Pegasi b 
has a mass of Mp = (0.46 ± 0.02) Mjup and an or- 
bital inclination between 79.6° and 82.2°, with the 
upper limit set by the non-detection of transits in 
photometric monitoring of the host star. By com- 
paring the cross-correlation signal obtained with 
a range of atmospheric models, we estimate mini- 
mum molecular abundances of VMR(CO) = 10"'* 
and VMR(H2 0) = SxlO"''. 

However, the absence of the signal in the third 
night of observations is puzzling. Based on a con- 
servative signal-to-noise estimate, there is a 20% 
probability of observing a signal less than 2.5a in 
the third night purely due to noise fiuctuations. 
Moreover, the stability of the CRIRES spectro- 
graph was much worse during the third night of 
observations, so that possible flat-field inaccura- 
cies affect a larger portion of the data. Further- 
more, residual cross-correlation signals from stel- 
lar features are present at the expected planet 
position in the night of October 25. We also 
investigated possible astrophysical causes of the 
lack of signal, i.e. the presence of a mild tem- 
perature inversion at the sub-stellar point or at- 
mospheric variability due to weather. However, 
none of the possible instrumental and astrophysi- 
cal causes can be confidentially linked to the non- 
detection, meaning that future observations of 51 
Pegasi will be essential in confirming the signal 
presented in this work. 

Ground-based high-resolution spectroscopy in 
the near-infrared is a powerful tool for investigat- 
ing hot Jupiter atmospheres. While the interpre- 
tation of low-resolution studies during transits or 
secondary eclipses is often controversial, due to the 
fact that molecular features superimpose across 
the same wavelength range, this technique unam- 
biguously identifies molecular signatures thanks to 
the use of a cross-correlation technique with tem- 
plate models. Moreover, it delivers the planet ra- 
dial velocities directly and does not require - when 



applied to dayside spectroscopy - a planet to tran- 
sit. The only requirement is that the planet must 
move significantly in radial velocity (by 10 - 20 
km s^^) during the observations. 

The main limitations of this technique are cur- 
rently the degeneracy between atmospheric struc- 
ture and molecular abundances, and the need for 
very bright targets. The former can in principle be 
solved by observing at multiple wavelengths and 
by detecting the signatures of different molecu- 
lar species, a goal that is within the reach of our 
current methods and of the available instrumen- 
tation. The latter will certainly benefit from the 
next-generation of telescopes with mirrors of 30 m 
or more, and possibly from spectrographs capable 
of observing larger wavelength ranges at once, as 
the total cross-correlation signal is approximately 
proportional to the square root of the number of 
strong molecular lines observed. 

Future developments of high-resolution spec- 
troscopy may involve the study of the shape of 
the cross-correlation peak during transit, capable 
of revealing atmospheric rotation and global cir- 
culation patterns. Furthermore, by accumulat- 
ing enough signal-to-noise, it will be possible to 
study the phase curve of these hot Jupiter and 
trace changes in the molecular abundances and in 
the atmospheric chemical processes between their 
morning and evening sides. Eventually, by em- 
ploying a 39-m class telescope such as the E-ELT, 
it may be possible to detect oxygen in the at- 
mosphere of an Earth-like planet o rbiting in the 



habitable zone of an M-dwarf star (Snellen et al 



20131 ). This would be a considerable jump forward 
towards the detection of biosignatures in Earth 
analogs. 
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Night 


Kp Vgys Phase Significance 

(km s^^) (km s^"'^) range 


2010-10-16 
2010-10-17 
2010-10-25 
16 -1- 17 


122.3+24;^ -25.7lio:3 0-36 0.42 4.5cr 
IMJtlll -33.2iio*3 0.60-0.66 4.7c7 

0.49-0.54 <2.5a 
134.1±1.8 -33.2±1.5 n.a. 5.9a 



Table 1: Details of the single-night detections of 51 Pegasi b. The planet orbital velocity, the systemic 
velocity, the observed range in planet orbital phase and the significance of the detection arc given. The last 
row shows the nights of October 16 and 17 combined. A global phase shift of Aip = 0.0095 is applied to the 
data in order to make the systemic velocity of the night of October 16 and 17 combined to coincide with the 
systemic velocity of 51 Pegasi. 



Molecule 


Kp 


Vsys 


Significance 




(km s~^) 


(km s"^) 




H2O 


136.9±2.8 


-34.9+-* 


4.1(7 


CO 


134.2±1.8 


-33.2±1.5 


5.7(7 



Table 2: Details of the cross-correlation signals obtained by employing models with pure water or pure carbon 
monoxide, for the nights of October 16 and 17 combined. The planet orbital velocity, the systemic velocity 
and the significance of the detection are given. 
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